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Abstract 

This paper concerns the derivation of a Fokker-Planck equation for the correlation of two high 
frequency wave fields propagating in two different random media. The mismatch between the 
random media need be small, on the order of the wavelength, and their correlation length need be 
large relative to the wavelength. The loss of correlation caused by the mismatch in the random 
media is quantified and the limit process for the phase difference is obtained. The derivation is 
based on a random Liouville equation to model high frequency correlations and on the method 
of characteristics to characterize mixing in the random Liouville equation. Applications of such 
correlation loss include the monitoring in time of random media and the analysis of time reversed 
waves in changing heterogeneous domains. 

1 Introduction 

The energy density of high frequency waves propagating in highly heterogeneous media can be 
modeled by a Fokker-Planck equation in the phase space, i.e., the space of positions and momenta. 
We refer to e.g. |2] for a mathematical derivation of such a macroscopic model for the wave energy 
density. The main assumption on the heterogeneous medium is that it is a random medium with a 
very large correlation length relative to the typical wavelength in the system. The Fokker-Planck 
equation may be seen as a highly-peaked-forward-scattering approximation to the radiative transfer 
equations, which are also used in the modeling of the energy density of waves in heterogeneous media 
when correlation length and wavelength are comparable; see e.g. [Tl [7| I13j. 

Such macroscopic models for waves in random media can more generally be used to quantify 
the correlation function of two wave fields propagating in possibly two different media 0. This has 
applications in the temporal monitoring of the statistical properties of random media as well as in 
the analysis of the refocusing of time reversed waves |3 01 E]- This paper analyzes the effect of 
changes in the random media on the two-field correlation. In the Fokker-Planck regime, the two-field 
correlation decays as the two media separate and we present a quantitative estimate of such a decay. 

Although the results in this paper generalize to fairly large classes of waves such as e.g., acoustic 
waves as in |2- electromagnetic waves, and elastic waves (see |13j ) we restrict ourselves to the 
case of a scalar Schrodinger equation to simplify. This models the effect of heterogeneities on a 
single particle represented by a quantum wave function. The changes in the heterogeneous medium 
considered here are sufficiently small so that the correlation of the two fields, one propagating in 
the unperturbed medium and the other one propagation in the perturbed medium, still satisfies a 
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random Liouville equation in the high frequency hmit. However, solutions of the Liouvihe equation 
are no longer real and have a complex phase. Their phase is responsible for the decoherence of the 
two wave fields, and is driven by the media mismatch. The evolution of the correlation function is 
described in terms of the particle position X^{t), momentum K^{t) and phase difference Z^{t). The 
techniques used in [2] are then applied to this new random Liouville equation to study the limit of 
the joint process [X^ {t) ^ [t] , {t)) . The limiting equation, of Fokker-Planck type, is obtained in 
the vanishing limit of the correlation length 5 of the heterogeneous medium. Its derivation is based 
on the mixing properties of the bi-characteristics of a random Hamiltonian, as in and on the 
mixing properties of a highly oscillatory functional of such bi-characteristics, which is the main new 
result obtained in this paper. In particular, we show that the phase difference Z^(t) converges to a 
Brownian motion on the real line. 

The rest of the paper is organized as follows. Section|21presents the random Schrodinger equations 
and the limiting random Liouville equation for the Wigner transform of the two wave fields, which 
is the Fourier transform in the difference variable of the correlation function of the two fields. We 
also review in Section 12.11 the results obtained in j^l while adapting them to the case of a random 
Schrodinger equation. Finally, in Section [2. 21 we derive the generalized Fokker-Planck equations (see 
1)2. 23() below) by formal analysis. Section |21 presents Theorem 13. 11 the main result of this paper that 
describes the limit of the joint process {X^{t), K\t), Z\t)). The Fokker-Planck equation is deduced 
from the limiting law (as the correlation length of the random medium goes to zero) . The proof of 
Theorem 13. II is presented in Section 0] 
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Fellowships. GB has been also supported by NSF grant DMS-0239097. 



2 The random momenta Liouville equation 

We consider the evolution of the correlation function of two solutions of the Schrodinger equation 
with a small mismatch between the random potentials but with the same initial data. The function 
satisfies 

ie^ + yA^, - Vs{x)A = (2.1) 

and the function cj)^ satisfies 

ie^ + -^4>e - mx) + eSs{x)](l)e = 0. (2.2) 

Both of the functions (j)e and "ipe satisfy initially 

V'.(0,a;) =0,(O,x) =(^§(x). (2.3) 

The family (j)^ is e-oscillatory and compact at infinity [S] . The random potentials V5 and 5"^ vary on 
a scale 6 that is much larger than the wave length e of the initial data but is much smaller than the 
overall propagation distance that is of the order 0(1): e<C(5<Cl. To keep a non-trivial correlation 
of ipe and (j)e the mismatch of the potentials has to be weak - hence the coefficient e in front of 
5"^. We will see that in order to produce an order one contribution we will have eventually to take 
Ss{x) = 5~^^'^S{x/5) making the overall strength of the mismatch be of the order 0(e/\/^). 
In order to study the correlation of tps and (ps we introduce the cross Wigner transform as 

We{t, x,k)=l e^>^-yA {t, x-'-l-)^,{t,x + f)j^. (2.4) 
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When (f)^ = ipi; then is real and is usually interpreted as the phase space energy density of the 
family (pei^)- However, the distribution W£(t,x,k) does not have to be real if (j)£ ^ ipe- Its phase 
measures the decoherence of the functions 4>£ and tpe. The basic properties of the Wigner transforms 
can be found in jBl . 

In order to obtain an equation for Wg we differentiate the Wigner transform with respect to time: 



dt 
1 

ie 



Jk-y 



Vs(x 



2 



V5[x + 



eSs (x + 



2 



X + 



ey\ dy 



(27r)" 



Passing to the limit e — > we obtain an equation for the distribution Ws{t,x,k), the weak limit of 
Wease^ 0: 

+k- V^Ws - VVsix) ■ VkWs = iS5{x)W5. 



dt 



(2.5) 



In order to obtain a non-trivial limit of Ws{t,x,k) as the correlation length 5 — > we choose the 
random potential Vs{x) and the mismatch Ss{x) to be of the form 



Then ()2.5() becomes 



£^,,.v.„..--Lw(|).v.,...-;=5(|)w. 

W{0,x,k) = Wo{x,k). 



(2.6) 



The initial data Wo{x,k) is simply the limit Wigner measure of the family cp^. Equation (|2.6|) is 
the starting point of our analysis. We note that if we consider the initial data for the Schrodinger 
equation as a mixture of states then the error bound for the approximation for the Wigner transform 
We by the solution of 1)2. 6|) may be estimated and the sequential limits e — > first, and 5^0 second 
may be replaced by a joint limit (e, 5) — > 0; see for details. This may be done in a certain 
region of the {e, 5)-plane that ensures that the scale separation e <C 5 <C 1 is kept under control. We 
will not pursue this avenue in this paper in order to avoid unnecessary technical complications. 



2.1 A review of the case in the absence of a mismatch 

We first recall the known results of (see also [S] for a longer time scale analysis) when there is 
no potential mismatch, that is, when 5 = and d > 3. We restrict our attention in this paper also 
to the case d > 3 though a generalization using the results of and jT|3] is possible. If 5" = then 
solution of ()2.6() is given explicitly in terms of the random characteristics {X^{t),K^{t)): define 

We have Ws{t,x,k) = Wo{X^ (t) , (t)) . The characteristic trajectories satisfy a limit theorem; the 
process K^{t) converges to a Brownian motion K(t) on the sphere {|A:| = |A;(0)|} and X^{t) converges 
in law to ^ 

X{t) =x+ [ K{s)ds. 
Jo 
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Before formulating the limit theorem we describe first the necessary assumptions on the random 
potentials V and S; we will not, of course, need the assumptions on S in the present section but we 
will need them later on and it is convenient to put them here. 

Let (r2,S,P) be a probability space, and let E denote the expectation with respect to P. Let 
us denote the joint process F{x,uj) = {V {x , uj) , S {x , u)) and assume that the random field F is 
stationary in the first variable. This means that for any shift x E M*^ and a collection of points 
xi, . . . ,Xn S the laws of {F{xi + x), . . . ,F{xn + x)) and (F(xi), . . . ,F(x„,)) are identical. In 
addition, we assume that E{S'(x)} = E{y(x)} = for all x G M'^, the realizations of V{x) and S{x) 
are P a.s. C^-smooth in x and they satisfy 

M := max ess-sup \d^ F{x,uj)\ < +00. (2.8) 

We suppose further that the random field F{x,uj) is strongly mixing in the uniform sense. More 
precisely, for any iZ > we let CJ^ and Cf^ be the cj-algebras generated by random variables F{x) for 
all X £ Mji and x G respectively. The uniform mixing coefficient between the cr-algebras is 

(P{p):=sup[\¥{B)-¥{B\A)\: R > 0, A e Cj,, B e 6%^^], (2.9) 

for all p > 0. We suppose that 0(p) decays faster than any power: for each p > 

hp := sup pP(t){p) < +CO. (2.10) 

The two-point spatial correlation 2x2 tensor of the random field F is denoted by R{y) and has 
components 

iZ^^(y)=E[y(0)l^(y)], ' (y) = E{V (0) S (y)} , i2^^(y) = E{5(0)F(y)}, R^' (y) = E{S{0)S{y)}. 
Note that (|2.in|) implies that for each p > 

hp := s^^P (1 + \y\y^'\d'^R{y)\ < +00. (2.11) 

We also assume that the correlation tensor R{y) is of the C°°-class and that 

R^^ (k) does not vanish identically on any hyperplane Hp = {k : {k ■ p) = 0}. (2.12) 

Here (k) = J R^^{x) exp(— i/c • x)dx is the power spectrum of V. 
Let us define the diffusion matrix Dmn by 

Drnn{k)=--j — — = / ^ — ^ "I, 71 = 1, . . . , d, k = k/\k\. 

(2.13) 

Then we have the following result. 

Theorem 2.1 |^ Let he the solution of ^2.6} with the initial data WQ{x,k) supported in a 
compact set away from k = 0; supp{Wq\ C S x A{M) with A{M) = {M~^ < \k\ < M} for some 
M > and a compact set S CM.'^. Let the function (j) satisfy 

|,,.V.^.£(z,,„.(.,i|) 

(j){0,x,k) = Wo{x,k). 
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Then, there exist two constants C{T) and > such that 



sup 

{t,x,k)&[0,T]xK 



mV^ {t,x,k)-(t){t,x,k) <C{T){l + \\Wo\\ci)5°° (2.15) 



for all compact sets K C A{M) = M'' x A{M). 
Note that 

^ ,,,r 1 f°° d^R^^(sk)^ , I d fdR^^{sk)\ , 

and thus the A'-process generated by (|2.14|) is indeed a diffusion process on a sphere k = const, 
or, equivalently, equations (|2.14|) for different values of \k\ are decoupled. It is easy to check that 
assumption (|2.12j) implies that the matrix D{k) has rank d — \ for each k G R'^jO}. It can be also 
shown that then equation ()2.14p is hypoelliptic on the manifold x {|A;| = k^} for each k^ > 0. 

2.2 A formal analysis of the momenta Liouville equation 

We first present a non-rigorous formal multiple scales analysis of ()2.6|) . which provides a short and 
relatively quick way to the correct limit. We introduce a multiple scales expansion 

Ws = W{t, X, k) + ^Wi{t, X, y, k) + 5W2{t, x,y,k) + ..., y = x/6 

and insert it into ()2.6|) . As usual we make an additional assumption that the leading order term 
W{t, X, k) is deterministic and does not depend on the fast scale variable y. In the leading order we 
obtain 

k ■ VyWi + eWi = VV (y) • VkW + iS (y) W. 

Here ^ > is an auxiliary regularizing parameter that we will send to zero at the end. Define the 
correctors Xj ^i^d rj as mean-zero solutions of 

dV 

k ■ Vyx, + ex, = ^ 

k-Vyri + 9r] = S{y). 

They are given explicitly by 

Jo ciyj 

and 

poo 

7]{y,k)= e^^'S{y- sk)ds. (2.17) 



JO 

The function Wi is given in terms of the correctors as 

d 



Wi{t, X, y,k) = J2 Xj{y, k)^^^^^j^ + ir?(y, k)W{t, x, k). 



The equation for W2 is 

dW 



dt 



+ k ■ V^W + k ■ VyW2 = VV{y) ■ VkWi + iS{y)Wi. 
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Averaging under the assumption that K{k ■ VyW2} = we obtain the fohowing closed equation for 
the leading order term W: 



dW 



+ k ■ V.VF = E {VV{y) ■ VkWi + iS{y)Wi} = Ji + Jn- 



(2.18) 



The two terms on the right side are computed using the explicit expressions (|2.16|) and (|2.17|) for 
the correctors. The first term may be split as 



with 



J/ = E {VV{y) ■ VkWi] = J] + J| 
dW{t,x,k)' 



Xmiy,k)- 



_d_ 

dki 



E 



dV 

oyj JO 



dkr. 



dV{y-sk)_^_,^^^\ dW(t,x,k) 



dyn 



dk„ 



_d_ 

dk^ 



dW{t,x,k) 
dkm. 



where the diffusion matrix Djm is given by (Eini). The term j| is 



Jf = E 



d 



dV , ^ d 



(y)^ k)Wit, X, k)]} =i 



dyj dkj 



d_ 

'dkj 



E 



dV 

9yj 



(y) I S{y - sk)e-'''dsjW{t,x,k) 



with the drift 



{E'^ik)Wit,x,k)) 



E^{k) 



dR^^isk) 



dxj 



ds. 



Now we look at the second term in the right side of (|2.18|) 

JH = ^{iS{y)Wi} = J]j + J]i 



(2.19) 



with 



Jh = E\iS{y)xm{y,k) 



dW{t,x,k) 



dkr. 



iE{S{y) 



dV{y-sk)_^_es^^\ dW{t,x,k) 



dyn 



dk„ 



,E,^^j^fWit,x,k) 



with 



Note that 



E" 



dkm 

°° dR^^jsk) 
dxi 



ds 



dR^^{-sk) 



dxj 



ds 



^ dR^^jsk) 

oo d-^j 



ds. 



J] + Jh = {E'^ik)W{t,x,k)) +zFm{k)^^^^^^ 



^(^- + + .(V, • E')W{t, X, k) = E,^-^^^^j^ + FW{t, X, k) 



with 



dkj 



OR^^jsk) 
dxj 



ds 



(2.20) 
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and 

F = Vk-E'= / sAR^^{sk)ds. (2.21) 
Jo 

The last term in (fTT!?]) is 

J]j = E{iS{z)r]{z,k)W{t,x,k)} = -K{k)W{t,x,k) 
with the absorption coefficient 

poo 

n{k) = / R^^{sk)ds. (2.22) 







Putting together all the terms above we get the equation for W: 

dW d ( dW\ 

+ k ■ = iE{k) ■ VkW + iF{k)W + -— Dmnik)^ - K{k)W. (2.23) 



Ot Okra V Okn 

If S and V are independent then F = E = and this simplifies to 

dW d ( dW\ 

^},.^.^W = -—\Drr.n{k)^\-K{k)W. (2.24) 



ot Okm V Clkn 

Unfortunately, the asymptotic expansion described in this section may not be justified rigorously. 
In the next section we present a rigorous derivation of the limit equation (|2.2H|) . 

3 The Liouville equation and the phase diffusion 
3.1 An example: the decorrelated case 

The purpose of this section is to obtain the limit equation 1)2.23(1 as the limit of l|2.6() . Recall that 
solution of (|2.6j) 

+k- V^Ws - -^VV (^) • VkWs = -^S ) Ws (3.1) 



dt ' ^ " Vs v<^y " " ^ V5 

may be obtained by the method of characteristics. Along the trajectories we have 

,X' 1 ^^^s(^) (3.2) 



dt ' dt y/s \ 5 J ' dt y/6 \ 6 

with the initial data 

X^{0)=x, K^{0) = k, Z\0) = 0. 
Then the solution of ()3.1() is given by 

W\t, X, k) = e'^'^'^Wo{X\t), K^t)). 

The results described in the previous section tell us that converges to a Brownian motion K{t) 
on the sphere, with the diffusion matrix Dmni^)-, and converges to its time integral: 



K^{t)^K{t), X\t) ^ X{t) = x+ [ K{s)ds. 

Jo 



We will show below that converges to a Brownian motion with the diffusion coefficient Hi{k) given 
by (|2.22() . In the simplest case when V and S are uncorrelated the Kolmogorov equation for the 
limit process {X{t), K{t), Z{t)) is 



dt ^ dkm V dkn J dz"^ ' 
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Recall that 

W\t, X, k) := E{W^{t, X, k)} = E,,fc,,=o {e'^'^''^Wo{X^{t), K^{t)) 
Therefore, in the limit 5^0, the function {t, x, k) converges to 

W{t, X, k) = g{t, x,k,z = 0). 
Here the function g satisfies the Kolmogorov equation 

g(0,x,k,z) =e"Wa{x,k). (3.4) 
It may be written as g{t, x, k, z) = e'^^q{t, x, k), where the function q satisfies 

q{0,x,k) =Wo{x,k). 

We see that actually W{t, x, k) = g{t, x, k, z = 0) = q{t, x, k). Note that ()3.5() is nothing but (|2.24p : 
this relates the approach of the present section to the formal result of the previous section. However, 
the Kolmogorov equation H3.3|) provides the description of the whole limit process Z(t) while H2.24() 
is just one of its reductions. 

3.2 The main result 

We will prove the following theorem. 

Theorem 3.1 The joint process {X^ (t), (t), (t)) converges in law in the limit (5 — > to the 
dijfusion process {X (t) , K (t) , Z (t)) with the joint generator 

£ct> = Dmn{k)^J-^ + [Dm{k) + D^{-k)\ + D{k)^ + E^{k)^ + E{k)^ -k-V,4> 

OkmOkn OkmOZ OZ"^ Okm OZ 

(3.6) 

with the coefficients 

Dmn{k) = - I R^l{sk))ds, D{k)= I R''''{sk)ds, (3.7) 







?VV ^„ n^.^ _ / uSS/ 

'"mn 





oo 

/poo 
R^isk)ds, Emik) = -J^ sARl^i 



sk)ds, (3.8) 



and 



POO 

E{k) = / sAR^^{sk)ds. (3.9) 
Jo 



The generator may be written slightly more compactly as 

i -a!: (-"(^Ol) .| («„(-.,£) .«,.)g-..v.,. ,3.10) 
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3.3 A formal computation of the limit 

We first present a formal computation, whicli is similar to the derivation in section 12^21 and leads to 
the generator ()3.6|) . We start with the Liouville equation including the phase variable 



|,..v.,--Lv.(|).v..--L.(|)|.„ 

and consider an asymptotic expansion of the form 

(j){t,x, k, z) = (i){t,x, z, k) + VS(l)i{t,x,y, z, k) + 502(i, x,y, z,k) + . . . , y = x/S. 
In the leading order we get: 



(3.11) 



k ■ VycPi + 001 = VV {y) ■Vk4> + S {y) 



dz ' 



As before > is a regularizing parameter that we will send to zero later. The leading order term 
(j){t, X, k) is assumed to be deterministic and independent of the fast variable y. Using the correctors 
Xj and r] given by 1)2. 16|) and 1)2. 17|) . respectively, we obtain an expression for cpi as 

(Pi[t,x,y,k,z) = Xj{y,k)— +r]{y,k) — . 

The equation for cf) is obtained in the same way as before, as a formal solvability condition for (j)2; it 
reads 

(3.12) 



-^ + k-V,cl) = EiVV{y)-Vk^i + S{y) 



The first term on the right is 



E{V{y)-VM = ^^^{vj{y) V^{y - sk)e-'' ds-^ + S{y - sk)e-'' ds^^ | 



d 



dk, [^"^'^^hkj^dk, 



d 



(3.13) 



in the limit 0^0. The second term in the right side of (|3.12|) is 



1 1/„ ^'<''=>'''a.„ 



+ D 



d 



dz"^ dz 

as ^ 0. Putting together 1)3. 13() and 1)3.14(1 we obtain 



R'J{-sk)ds^+D 



dz"^ 



(3.14) 



This is nothing but the Kolmogorov equation for the process with the generator (|3.1U() . 
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4 Proof of Theorem IM.ll 



Theorem 13. II is a simple corollary of the following proposition. Let us first introduce some notation. 
Given a function G e C^([0, +00) x M^'^) and t > let us introduce 



Nt{G) = Git, Xit),Kit), Z{t)) - G{0, XiO), K{0), Z(0)) - {8, + C)G{p, X{p),K{p), Z{p))dp 



with the operator C defined in We also denote by Wl := \ {0} and M^"' := M"' x to 

avoid the singular point k = 0. We let C := C([0, +00); M*^ x x M) be the set of continuous 
paths of {X{t),K{t), Z{t)). For any u < v denote by Ai^ the cr-algebra of subsets of C generated by 
{X{t),K{t),Z{t)), t £ [u,v]. We write M" := Mq and M for the a algebra of Borel subsets of C. 
It coincides with the smallest cr-algebra that contains all t > 0. We define C{T,M) as the set 
of paths vr G C so that both {2M)-^ < \K{t)\ < 2M, and 



i 

X{t) - X{u) + j K{s)ds = 0, for all < u < t < T. 



We denote by ^ ^ the probability measure on C(T, M) induced by the trajectories of (|3.2|) and by 



the corresponding expectation. 



Proposition 4.1 Suppose that {x,k) G A{M) =R'^ x {M"^ < |A;| < M} for some M > and that 
a test function G Cft((M^'^)"') is non-negative. Let 70 G (0, 1/2) and let < ti < • • • < tn < 
t < u < T. Assume further that u — t > 6^^° . Then, there exist constants 71 > 0, C(T) such that for 
any function G G C^{[T^,T] x M^rf x R) we have 



e; 



{ [iV4G) - Nt{G)] C} I < C{T)6^- {u - t)\\G\U (e^,,,/ + 



(4.1) 



Here C(vr) := C{X{ti), K{ti), Z{ti), . . . , X{tn), K{tn), Z{tn)), and tt = {X (t) , K (t) , Z (t)) is any con- 
tinuous path. The choice of the constants 71, C does not depend on {x,k,z), 5 G (0, 1], the times 
ti, . . . ,tn, 7*, T, u, t, or the test function G. 

Proof of Theorem 13. II Theorem 13. H is a simple consequence of Proposition 14.11 Let (/>o(a;, k, z) 
be a test function and let the function (f)[t, x, k, z) solve the initial value problem 



(4.2) 



(j){0, X, k, z) = (j)o{x, k, z). 

We apply Proposition 14.11 with G{t, x, k, z) = (j){u — t,x,k,z), t = 5^ and u > 6'^ with 1/2 < 7 < 1 
and take C = 1- It follows from 1)4.1(1 that 

E^,^, [MXiu),Kiu), Z{u)) - 4>{u - 6\ X{5^),K{5^), Z(5^))] | < C||G||4<5^^ (4.3) 

Using the fact that is a smooth function and l/2<7<lwe conclude that 

Wi^f,A<i)Q{X\u),K^{u),Z^{u)) -4){u,x,k,z)\ < C||G||4(5^i. (4.4) 

The conclusion of Theorem 13. II now follows. □ 
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4.1 The proof of Proposition (4.11 

The proof of Proposition 14.11 is technical and uses the ideas of |2llSliniEl- However, the present 
situation is much simpler than in the aforementioned papers as we already know from Theorem 12. II 
that the process {X^,K^) converges to a process with the generator 

This means that the law of the process will be close to that of the limit X{t); and in particular 
X^{t) does not approach a narrow tube around its past trajectory with a probability very close to 
one. This potential return was a major obstacle in the proofs in pi 151 1^ ITU]. 

The strategy of the proof is as follows. We will need to deal with objects of, say, the form 

E \GiX^{s),K\s), Z\s))S ( ^\ V ( ^\ \ (4.5) 



with s < s' but the difference s' — s small. Then we will consider a slight pullback in time 

a{s) = s - 5^-^' , (4.6) 
with a sufficiently small 7 > and a linearization 

L{a, s) = X^{a) - (s - a)K^{a). (4.7) 
The characteristic equations (j3.2j) allow us to estimate the difference between (|4.5|) and 

E {G{xHa),KHa), Z'{a))S (^) V } , (4.8) 

and show that it is small. However, the latter expectation approximately splits: 

^\G{X\alK\alZ\a))S (^) V (^) } - (4-9) 



E \G{X\a), K\a), Z\a))^ R^^ 



L{a, s') — L{a, s) 



The reason for the expectation splitting is that the argument of the function G in 1)4. 8|) depends only 
on the potential in a tube close to the trajectory X^{t) until the the time t = a, while L{a, s) and 
L((T, s') are at distance much larger than 6 from this tube with a probability close to one. Hence, the 

values of G{X^{a), K^{a), Z^{a)) and, say, of V ( — ^ — \ are nearly independent and expectation 



<5 

1)4. 8() splits. This is formalized by the following mixing lemma. 

For any t > we denote by !Ft the u-algebra generated by {X^{s), K^{s), Z^{s)), s < t. Here we 
suppress, for the sake of abbreviation, writing the initial data in the notation of the trajectory. We 
assume that Xi,X2 ■ (M x M'^ x M'^^)^ — > M are certain continuous functions, Q is a random variable 
and 5i, 52 are M'^-valued random vectors. We suppose further that Q, 51, (?2; are .Fj-measurable, while 
Xi , X2 are random fields of the form 

X^{x)=Xi {F{x),V^F{x),VlF{x)) , 

where, as before we denote for brevity F = {V, S). We also let 

?7(0i, ^2) := E [Xi(0i)X2(02)] , 01,62 eR'^ (4.10) 

and recall that (j){r) is the mixing coefficient defined in 1)2. 9|1 . The following mixing lemma from 
is formalizing the expectation splitting. 
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Lemma 4.2 (i) Assume that r,t >0 and 



inf 

u<t 



r 



(4.11) 



-a.s. on the set Q ^ for i = 1,2. Then we have 



E 



Xi{gi)X2{92)Q -K[U{g^,g2)Q] 



< 2</. ( -j \\Xi\\L^\\X2\\L^\\Q\\Liin)- (4.12) 



(ii) Let EXi(O) = and assume that §2 satisfies 



inf 

u<t 



91 



X\u) 



6 



> 



r + ri 



(4.13) 



and \gi — (72I > fi^ ^ for some ri > 0, F-a.s. on the event Q ^ 0. Then, we have 



E 



26 



Xi{gi)X2{g2)Q -nU{gug2)Q] <C0i/2(^)0i/2(ri) ||Xi||i^||X2||L^||Q||ii(m (4.14) 
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for some absolute constant C > 0. Here the function U is given by \4^.1U^ . 

We proceed now with the proof of Proposition 14. II Let G{k, z) be a sufficiently smooth function. 
We will establish the approximate martingale property 1)4. 1() for G. It suffices to consider functions 
of the form G{z, k) = g{z)r{k). The characteristic equations 



and 



imply that we have 



G{K'{u),Z'{u))-G{K\t),Z'{t)) 



(4.15) 



1 



g'{zHs))S 



X\s) 
S 



r{K\s))+g{Z'{s))r,{K\s))V, 



X'{s) 
6 



ds. 



Here and below, we use the notation Vj for Vx^ and similarly 5"^ for 5^^. to simplify. In order to be 
able to make a backward step a as in 1)4. 6|) we split the above integral as 



G{K\u), Z^{u)) - G{K^{t), Z\t)) =A + B 



(4.16) 



with 



A 



g'{z'{s))r{K\s))S 



xHs] 



+ g{Z\s))r,{K\s))V, 



X'{s) 



ds 



and 



B 



g'{Z\s))r{K'is))S 



X'(s] 



+ g{Z'{s))r,{K\s))V, 



X\s] 



ds = Bi + B2. 



The first term is small: 



\A\ < C7(^i/2-7||G<||. 
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provided that 7 < 1/2. The term B will be analyzed with the aforementioned ideas of a backward 
step and linearization, and using the mixing lemma. The terms Bi and B2 may be written as 

i?i = ^ r {g'{Z'{a{s)))r{K\a{s)) + [g' {Z\s))r{K' (s)) - g' {Z\a{s)))r{K' {a{s))]} 
xs(^^^]ds = I + II (4.17) 



and 



B2 = ^ r {g{Z'{a{s)))r,{K'{a{s)) + [g{Z' {s))r^{K' (s)) - g{Z\a{s)))r,{K\a{s))]} 

(4.18) 



X 1A- 



ds = III + IV 



with 



and 



while 



and 



V6 



g\z\a{s)))r[K%a[s)))S 



ds 



1 



II = / [g'{Z\s))r{K\s)) - g' {Z' {a{s)))r{K\a{s))]S 



III 



1 



1 

71 



g{Z\a{s)))r,{K'{a{s))V, 



X\s) 



ds 



V6 



IV = ^ _ [g{Z\s))r,iK'is))-g{Z\a{s)))r,{K\ais))]V, 



6 



(4.19) 



ds (4.20) 



(4.21) 



ds. (4.22) 



We summarize the contributions of each of the terms above in the following lemma. 

Lemma 4.3 There exists a constant a > so that the terms I, II, III and IV satisfy the following 



estimates: 



E 




OAR'^^ (eK\s)) deg'{Z'{s))r{K'{s))ds 



t JO 

<C6''\\G\U [||CI|oo + ie{c}] (u-t), 

U fOO 



c 



E 



II 




g"{Z\s))r{K%s))R''^eK'is))d6ds 



(4.23) 



(4.24) 



t JO 




t JO 



g'{Z\s))rj{K\s))Rf {eK\s))deds \ ( 

ru /"oo 



nCXD 
eg{Z\s))rj{K\s))ARj^ (oK^s)^ dOds 

<c5-\\Gh(\\aoo+m) {u-t) 



<C'^"l|G'||3 [llCI|oo + i^{c}] {u-t), 

(4.25) 



E{iyc} 



g'{Z\s))rj{K\s))Rf [-eK\s)\ dOds 



sv 



(4.26) 



t 



+ 



u 00 

J J giZ'{s))rjm{K'{s))RlJ (OK^s)) dOdsC 



t 



<C5^Ghmoo + E{Q]{u-t). 
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It remains now only to prove Lemma 14.31 as the four individual contributions above combine to the 
operator L in (|3.6() . 



4.2 The estimate for / 

We first recall the linear approximation l|4.7j) and define the interpolation 

R{v, a, s) = (1 - v)L{a, s) + vX\s). 
This allows us to write a linear approximation for S as 

_ ^ [ R{l,a,s)\ _ ^ [ R{<d,a,s)\ ^ 1 ^ [ R{v,a, s) ^ ^ 



S ' =S ^ ' ' ^ +7 / {Xt{s)-U{a,s))dv 



S ( ^) +U'sJ ] iXtis) - Ha, s))dv. 



5 J \ ^ J \ 6 J 6 

1 
'6 

Now we split I as 

/ = Ji + J2 (4.27) 



according to the above, with 

1 

T^Jt 



Ji = ^ r 9'{Z'{a{s)))r{K'{a{s))S (^^^ ) ds (4.28) 



and 



The term Ji is ready for an application of the mixing lemma: the arguments of the function G (that 
is, of g and r), and of the field S are separated by a distance of the order 0{6^~'^) that is much 
larger than 5, with a probability close to one. In order to make this statement precise we introduce 
a stopping time that ensures that until ts the trajectory X^{t) "goes forward" and does not come 
back to its past. 

Let < ei < £2 < 1/2, £3 G (0, 1/2 — £2), £4 E (1/2, 1 — £1 — £2) be small positive constants and 

set 

AT = [S-'^], p = [6-'% q = p [5-^=*], A^i = Np [6-^^]. (4.30) 

The requirement is that Ei, i € {1,2,3} should be sufficiently small and £4 is bigger than 1/2, less 
than one and can be made as close to one as we would need it. It is important that £1 < £2 so that 
N ^ p when 6 <^ I. We introduce the following {M*)t>o stopping times. Let t^^ := kp~^ be a 
mesh of times, and vr G C be a path. We define the "violent turn" stopping time 



t > : for some A; > we have t E 



t'ilt^l) and (4.31) 



Vsiir) := inf 

i^(4-2,) . m <i-^,ork (4^) - i-) . m < 1 - i 

where by convention we set K{—l/p) := K(0). Note that with the above choice of £4 we have 
k (^t'^^ - l/iVi) • -^(4^^) > l-l/N, provided that 6 G {0,6o] and 60 is sufficiently small. The 
stopping time Vs is triggered when the trajectory performs a sudden turn; this is undesirable as the 
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trajectory may then return back to the region it has akeady visited and create correlations with the 
past. 

For each t > 0, we denote by Xtiir) := [J X (s;7r) the trace of the spatial component of the 



0<s<t 



path TT up to time t, and by Xt((?;vr) := [x : dist (x, j£t(7r)) < 1/q] a tubular region around the path. 
We introduce the stopping time 



UsiTT) := inf 



t > : 3 A; > 1 and t G 45i) ^o^' ^^ich X{t) G {q) 



(4.32) 



It is associated with the return of the X component of the trajectory to the tube around its past; 
this is again an undesirable way to create correlations with the past. Finally, we set the stopping 
time 

Tsiir) := Vsiir) A Us{Tr). (4.33) 



Lemma 4.4 The probability of the event [rg < T] for a fixed T > goes to zero, as 5 — > 0.' 
there exists qq > so that 

^{[rs <T]} <C{T)6''°. (4.34) 



We apply part (i) of Lemma 14.21 to E{Ji(^} with 

Xi{x) = S{x), X2 = l, Q = g'{Z^{a))r{K^{a{s)))l[Ts>T]C, gi 



L{a,s) 



Note that gi and Q are both measurable. It follows from the definition of the stopping time ts 
that when Q ^ then the linearization also stays away from the past trajectory: for all < p < a 
we have 



Lia,s)-X\p) 



(4.35) 



and hence 



inf 

0<p<o-(s) 



91 



r 

> - 
- 6 



with r = C6^ ^. We also decompose Ji according to whether the stopping time has occurred or not 
before time T: 



1 



g'{Z'{a{s)))r{K'{a{s)))l[Ts>T]S 



L(cr, s) 



6 



ds 



+ -= / g'{Z\a{s)))r{K\a{smi - 1[ts > T])S 



ds = Jii + Ji2. 



However, ()4.34l) implies that 



E 



{\ji2C\} <C5'"^\\G\\2\\C\\oo{u-t) 



(4.36) 



so we have to deal only with Jn. Using the mixing lemma as above, with the point separation as in 
(jO^ . and the fact that E[S'(x)] = (whence [/ = in ^J^) we estimate 



E JnC 



< {C6-^) (u - t)||5||i ||r||oE[C] < CSWGhnCKu - t). 



(4.37) 



15 



Here 0(r) is the mixing coefficient that decays faster than any power of r; see (|2.1()j) . We conclude 
that 



E JiC 



< C5\\GhE[C]{u -t)\\+ Cd'^'^WGUWooiu - t) 



(4.38) 



so the term Ji produces only a small contribution. 

Now we estimate the second term J2 in (|4.27j) : it is given explicitly by (|4.29j) . We split it further 
by using the next order expansion 



5,; 



R{v, a, s) 
S 



S^ 



L(a,s)\ 1 



Sij 



This leads to the corresponding expression J2 = J21 + J22 with 



J- 



1 



21 



53/2 



g'{Z\a{s)))riK'iais)))S. 



L{a,s) 



{Xf{s)- L,{a,s))ds 



and 



U 1 V 



J: 



22 



J5/2 



g'{Z\a{s)))riK'iais)))S., 



R{e,a,s] 



6 



t+5'^—i 



X (Xfis) - Li{a,s)){XUs) - Lj{a,s))dedvds. 



Note that the characteristic equations and the definition 1)4. 6|) of the time a{s) imply that 

\L{a, s) - X^{s)\ < ^52-27-1/2 ^ (^7^3/2-27^ (4 39) 

It follows that 



E 



J22C}| < C||G||2E{C}<5-5/253-47(^ _ < C 5^^^-^^ WGhnO {u - t), (4.40) 



which is small if 7 < 1/8. This is an important characteristic feature of the weak coupling limit; 
after several linearizations the remainder becomes deterministically small while the linearized terms 
may be controlled with the mixing lemma. 

Next, we look at J21 that is the only potentially surviving (not small) in the limit 5 — > 
contribution to /: to do so we write, using the evolution equation for X^ and a further linearization 
for the function Vi: 



Xf{s)-LM,s) = - f\Kf{v)-Kf{a)]dv = -^ f f V, 

Ja yd Ja JcF 

'X\p) 



1 

- f 



{s - p)V, 



1 



53/2 



(s-p) 

s rl 



V, 



6 



dp 



Vi 



1 

~7sJa 

L{a, p) 



is - P)V 
dp = — 



5 

L{cr, p) 
6 

1 r 

71 J a 



dpdv 
dp 



is - p)V, 



L{a, p) 



dp 



a JO 



(s - p)Vir, 



Riv,a, p) 



Xfnip) - Lmicr,p)^ dvdp. 



(4.41) 



This produces a further split 



J21 — •^21 + "^21 
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with 



A = f\s - p)g'{Z\a{s)))r{K\a{s)))S. V, dpds 



and 

u si 



im 



J2i = -p J J j{s-p)g'{Z\<y{s)))r{K\a{s)))S,[^^]V. 

t+5^-~l cr 
X i^X^rnip) - Lrn{(T,p)^ dvdpds. 

Note that we have hnearized enough to achieve a deterministic estimate 

\A\ < C6-^S^-^6^-^6^/^-^\\G\\i{u -t) = CS^/^-^^\\G\\i{u - t) (4.42) 

which is small if 7 < 1/6. Hence, J21 remains the only potentially contributing term in I: it is 
analyzed with the help of the mixing lemma. We take 

L{a,s) L{a,p) 

91 = — ^ — , 92 = — ^ — , Xi{x)=Si{x), X2[x)=Vi{x), r = p-a, n = s - p, 

and 

Q = g'{Z'{a))r{K'{a{s)))l[Ts>T]C. 

Let us denote 



Observe that 



Rf''{x) = -E{Si{y)Vj{x + y)}. 



f/ix) = -E{S,iy)V,ix + y)} = -^E{S,iy)V{x + y)} = -^E {S{y)V,{x + y)} = 



E 



The mixing lemma implies that 

u s 

4iC] + ^ / I{s-P)E l^g'{zHa{s)))r{K'{a{s))) [-R^ j^ Ha, p) - ^a, s) ^^^^ | ^^^^ 



since the mixing coefficient (f) is rapidly decaying. It follows that 



E 



4iC] - ^ f\s - P)E [g'{Z'{a{s)))r{K\a{s))) (a/?^^ (^^4^)) } ^^^^ 



<C5\\G\\iM0{u-t) 
The integral above from cr to s may be massaged as 
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Observe that 



\K\s)-K\a)\ < Cb^l'^-^ 



and the function R^^ is smooth and rapidly decaying. This allows us to replace a in the argument 
of by s. The same can be done with the functions g and r; we conclude that 



E 



7I 



t JO 



9AR^^ (9K'{s)) deg'{Z'{s))r{K'{s))ds 



c 



< (:75^i||G||4e|c} (u-t) (4.43) 



with some 71 > 0. Therefore, putting all the work of this section together, see (|4.:-{8j) . (|4.4()|) and 
1)4. 42() . we also have 



E 



/- / / Ql^R'"" [QK\s)]dBg\z\s))r{K\s))ds 
t Jo 



c 



This proves the estimate (|4.23|) in Lemma 14.31 



|oo+e{c}] (u-t). 

(4.44) 



4.3 Estimate for // 

We now look at the term // given by H4.2(J() and split it further as 

/I = i= r [g'iZ'{s))r{K'is)) - g' {Z' {a{s)))riK' {ais))]S ] ds 



1 

5 



i+<5l-T cr(s) 
u s 



+ 



6 



g"{Z\p))riK\p))S ( ^) S ) dpds 

g'{Z\p))r,{K\p))V, ( ^] S ] dpds, 



The estimation of // is very similar to that of / both in spirit and in mechanics but is even somewhat 
simpler since as all we have to justify is the replacement of the arguments of S and Vj by the 
corresponding linear approximation from the time (t(s). This is done as in the previous section with 
the help of the mixing lemma and linearization and leads to 



e{//c}-e{//'c}| <C5^1|G'||3 [||C||oo + ]E{C}] {u-t) 



(4.45) 



with 



g"{z\a))r{K\a))S 



-5 1 



L{a, p) 



S 



L{a,s) 



u s 



s 



g'{Z\a))r,{K\<y))V, 
The mixing lemma and rapid decay of the mixing coefficient imply that 



dpds 



dpds = III +Ih- 







e 


< 




_ V 



g"{Z\a))r{K\a))R'^ 



ssfL{a,s) - L{a, p) 



dpds > C, 



< C5\\GhE{C}{u-t) 
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The inner integral may be re-written as 



g"{Z'{a))r{K\a))R' 



,Ss(L{(j,s) -L{a,p) 



dp= lg"[Z\a))r{K\a))R 



ss 



{p-s)K'{a) 



dp 



<5-T 



g"{z\a))r{K\cy))R^'{eK\a))de 



Using the rapid decay of R^^ , smoothness of G and closeness of s and a we obtain 



E 



Ih 



t Jo 

Similarly, we have for II2 



g"{Z'{s))r{K\s))R^^{eK\s))deds \ ( 



<C6^'\\GhE{c}in-t). (4.46) 



E 




The inner integral may be re-written as 



< C5\\GhE{Q{u-t) 



g'iZ\a))r,iK\a))R) 



vs 



L{a,s) - L{a, p) 



dp 



g'{Z'{a))r,{K\a))R]^{-9K\a))de 



VSi 



5—< 



g'{Z\<y))rj{K\a))Rf{eK\a))dB 



Using the rapid decay of i?"^"^, smoothness of G and closeness of s and a we obtain 



E 



Ih 




t JO 



g'{Z\s))r,{K\s))R^''{eK'{s))deds ( 



Putting (|4.45j) . (|4.46j) and (|4.47jl together we obtain (|4.24|) . 
4.4 Estimate of /// 

We look at the third term in (|4.18|) given by (|4.21j) 

III = ^ /" g{Z'{a{s)))r,{K'{a{s))V, 



<G5^'\\Gh \\C\\^ + E\C\ (u-t). 

(4.47) 



ds. 



(4.48) 



This term is similar to / in that before using the mixing lemma we have to expand a little bit: 



L{a,s) 



<5 J 

Accordingly we split III as /// = I Hi + III2 with 



+ 7/ V,J^^^)iXlis)-L,ia,s))dv. 



nil 



1 



g{Z\a{s)))rAK\a{s)m 



L{a,s) 



ds 



(4.49) 
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and 



Ilh 



p/2 



j+5i~7 Jo 



g{Z\a{s)))r,{K\a{s))V,k 



(X^is) - Lk{a, s))dvds. (4.50) 



The expectation of the first term on the event [ts > T] is small by the mixing lemma because the 
points X^{a) and L{a,s) are at distance of order 5^~'^ , 'E[Vj{x)] = 0, and the mixing coefficient is 
rapidly decaying. On the other hand, P[r5 < T] is small according to Lemma 14.41 We conclude that 

E {lIhC}\ < Cd\\Gh (llClloo + E[C]) {u - t). (4.51) 
In order to estimate III2 we write 



R{v, a, s) 



5 5 



Vjkm 



R{9,a,s] 
6 



iXiis)-Lmic7,s))de, 



which splits III2 = III21 + III22 (note that the variable integrates out in //I21): 



HI'. 



1 



21 



g{Z'{a{s)))r,{K'{a{s))V,k 



L{a,s) 



and 



22 



^5/2 



- . - 



iXiis)-Lkia,s))ds, 



R{e,a,s] 



g{Z'{a{s)))r,{K\a{s))V,kr. 
X (Xi(s) - Lm{a, s)){Xi{s) - Lfc(a, s))dedvds. 
We have linearized sufficiently to make III22 be deterministically small because of (|4..S9)) : 

\in22\<C5-'''H^~^^\\G\\2{u-t). 

This leaves us with III21 to take care of. This we do with the help of 1)4. 41() 

'x'ipY 



(4.52) 



(4.53) 



(4.54) 



Xf{s)-ma,s) = - r[K{v)-K{a)]dv = -^ f f V, 



6 



dpdv 



{s - pm 



dp 
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X'ip) 



6 



1 

~7sJa 

L{a, p) 



is - p)Vj 



L{a, p) 



dp 



S 



1 



S fl 




is - p)Vj 



a JO 



R{v,a, p) 



X^ (p) - Lrn{cF, p) ] dvdp. 



that allows us to decompose III21 = III21 + ^-^-^21 with 



and 



21 



ml 



1 

"52 



t+(51-T Jc7 



[s - p)g{Z\a{s)))r,{K\a{s))V,k 



L(cr, s) 



1 



ff 

t+i^—y J a JO 



{s- p)g{Z\a{s)))r,{K'{a{s))V,k 
X {x^{p) - Lm{(T,p)^ dvdpds. 



V: 



L{a,s) 



L{a, p) 



Vjm 



dpds, (4.55) 



R{v, a, p) 



(4.56) 
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Again, III21 is deterministically small because of (|4.lffl)) : 

< CS-^6^-^^6^/^-^\\G\\2{u -t)< C6^/^-^^\\G\\2{u - t). 
Now, for III21 we first compute 



mk{x + y)Vjiy)} 



92 



dxjdxj^ 

Q3 



E{Vix + y)V,iy)} 



dxjdx, 



-E{Vj{x + y)V{y)} 



E{V{x + y)V{y)} = -AR]^^{x) 



dx'^dxk 



and use the mixing lemma to conclude that 

|e 

< C6\\G\\iE[C]{u-t). 



As before we change variables in the inner integrals above and replace the argument a{s) of smooth 
functions appearing above by s, as well replacing the limits of integration by their values in the limit 
5 — > 0, to conclude that 



E , 

't Jo 

< C(^||G||4E[C](n-t). 



eg{Z'{s))rj{K\s))ARj^ [-eK'{s)\ ddds 



Overall, the work of this section implies that 



E 



'n poo 

+ 1 J eg{Z\s))rj{K^{s))AR]^ [eK\s)]d9ds 



(4.57) 



<c6^G\u[\\c\\oo+m]{n-t) 

which is nothing but ()4.25|) . 



4.5 Estimate of IV 

We are now down to the last term IV in H4.18() that is estimated as // with the help of the mixing 
lemma and no additional expansions: 



u s 



IV 



1 



g'{Z\p))r,{KHp))S ( ^) V, (^] dpds 



1 



u s 



g{Z'ip))r,m{K'ip))Vm ( V, ( ] dpds. 



6 



f+<5l-T O-(s) 



First, we observe using the mixing lemma and smoothness of G that (compare to (|4.4 



E 



(4.58) 



IVc}-E[lV'c}\<G6^-\\Gh [||C||oo + E{C}] (u-t) (4.59) 
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with 



g'{Z'{a))r,{K\a))S 



L{a,p)\ ( L{a,s 



dpds 



(4.60) 



u s 



+ 



1 



g{Z\a))rjm{K\a))Vr, 



dpds = IV1+IV2. 



Now, we have by the mixing lemma 

u s 



g\Z'{a))rj{K\a))R) 



L{a,s) - L{a,p) 



dpdsQ 



<C5'^\\Ghm\oo + E{Q]{u-t) 



and thus 



g'{Z\s))rj{K\s))Rf [-eK\s)] dOdsQ 



sv 

'j 



t 



< C(5°||G||2 



+ E{C}){u-t). 

(4.61) 



Similarly, for the other contribution we have 



giZ\a))rj^iK\a))R 



mj 



L{a,s) - L{a,p) 



dpds( 



<C5-||G||2f||C||oo + ^{C}) {u-t) 



and 



ETOC} + J j g{Z\s))rj^{K\s))R^ [9K\s)]dedsC 
t 



<C6^G\\2(\\C\\oo + E{C}){u-t 



(4.62) 

Together, ()4.6ip and ()4.62p imply (|4.26p . This finishes the proof of Lemma 14.31 and thus of Proposi- 
tion as well. □ 
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